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Abstract The area of stable motion for fictitious Trojan asteroids around Uranus' equila- 
teral equilibrium points is investigated with respect to the inclination of the asteroid's orbit 
to determine the size of the regions and their shape. For this task we used the results of 
extensive numerical integrations of orbits for a grid of initial conditions around the points 
L4 and L5, and analyzed the stability of the individual orbits. Our basic dynamical model was 
the Outer Solar System (Jupiter, Saturn, Uranus and Neptune). We integrated the equations 
of motion of fictitious Trojans in the vicinity of the stable equilibrium points for selected 
orbits up to the age of the Solar system of 5 x 10 9 years. One experiment has been undertaken 
for cuts through the Lagrange points for fixed values of the inclinations, while the semimajor 
axes were varied. The extension of the stable region with respect to the initial semimajor axis 
lies between 19.05 < a < 19.3 AU but depends on the initial inclination. In another run the 
inclination of the asteroids' orbit was varied in the range 0° < i < 60° and the semimajor 
axes were fixed. It turned out that only four 'windows' of stable orbits survive: these are the 
orbits for the initial inclinations 0° < i < T, 9° < i < 13°, 31° < i < 36° and 38° < i < 50°. 
We postulate the existence of at least some Trojans around the Uranus Lagrange points for 
the stability window at small and also high inclinations. 

Keywords Uranus • Trojan asteroids ■ stability area 



1 Introduction 

The first discovery of a Jupiter Trojan in 1906 (Achilles by Max Wolf in Heidelberg) proved 
that the equilateral equilibrium points in a simplified dynamical model Sun - planet - mass- 
less body are not only of hypothetical interest. Ever since many of such Trojan asteroids of 
Jupiter have been found and now we have knowledge of several thousands of objects in the 
1 : 1 mean motion resonance with Jupiter. Several investigations show the symmetry of these 
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two st a ble equilibri um points not only in the restricted three body problem dSandor et alj 
d2002h : lErdil dl988h ) but also in the realistic dynamical model of th e Outer Solar System 
(OSS) consisting of Jupit e r, Saturn, Uranus and Neptune, like in e .g. [Dvorak and Sc hwar: 
d2005h: iFreistetted d2006h : ISchwarz et all d2004) ; iNesvornvl d2000h : iRobutel et all ( 1200: ' 
ITsiganis et al.l ( 12005; 



But we have also evidence for the existence of asteroids in orbit around the equilateral 
equilibr ium points of Neptu ne and also of Mars. Since the discovery of the first Neptune 
Trojan dChiang et all d2003h ) — the minor planet 2001 QR322 around its equilibrium point 
Lg - a lot of work has been dedicated to understand the dynamics of Neptune Trojans. On one 
hand the stability of hypothetical bodies was studied in the model of the Outer Solar System, 
on the other hand long term investigations tried to understand how Trojans can survive, or 
even get captured in these kind of orbits when one allows the large planets to migrate. 
An extensive stability study of the Trojans for the four gas giants has been undertaken by 
INesvornv and Donea d2002h in different dynamical models. The results of the realistic n- 
body simulations for clouds of hypothetical Trojans around the Lagrange point L4 showed 
how fast the depletion of asteroids is acting for Saturn, Uranus and Neptune. For several 
hundreds of bodies, different eccentricities and inclinations as well as semimajor axes the 
integrations were carried out over 4 Gyrs. Saturn Trojans turned out to be unstable already 
after several 10 5 years, Uranus Trojans were also found to be unstable but only after several 
million years, whereas Neptune's Trojans mostly survived for the whole integra tion time of 
4 Gyrs. In a similar study for Uranus and Neptune Trojans [Marzari et alj d2003h determined 
the diffusion speed of Trojans and derived even expressions for the secular frequencies using 
their numerical results. Both studies were undertaken for initial inclinations of the Trojans 
up to i = 30° and showed quite similar results for the dependence on the inclination. 



Taking into account the early evolution of the Solar System the migratio n processes need 
to be included. This has been done in a paper by iKortenkamp et al.l ( 120041) for the Neptune 
Trojans where it has been shown that only 5% of an initial Trojan population could survive 
this early stage when Jupiter migrated inward (5.4 to 5.2 AU), and Saturn (8.7 to 9.5 AU), 
Uranus (16.2 to 19.2 AU) and Neptune (23 to 29 AU) migrated out wards. Three different 
theorie s for the origin of Neptune Trojans are developed and tested by Chiang an d Lithwickl 
d2005h . they use the results to constrain the circ umstances of the formation of plan ets and 
the accretion of T rojans. In a more recen t study INesvornv and Vokrouhlickvl d2009h in the 
NICE model (e.g. ITsiganis et al.1 (2005b)) - where Uranus and Neptune could interchange 
their location in the Solar System and Jupiter and Saturn go through the 2:1 Mean Motion 
Resonance - it was shown that the Neptune T rojans were captured by a process quite sim- 
ilar to the chaotic capture of Jupiter Trojans dMorbidelli et al.l ((2005)). A similar process 
can be assumed for the Uranus Trojans, but up to now no Uranus Trojans have been de- 
tected which can be ex plaine d by the results of the p ure dynamical studies mentioned before 
jNesvornv and Done 2 d2002h :lMaT zari et al.1 J2003)). These investigations stop ped at an in- 
clinati on of i = 30° ; the discovery of a highly-inclined Neptune Trojan ( Shep pard and Truiillol 
d2006h ) raised the question of the gener al stability of Troja ns at high inclinations. In a recent 
study concerning the Neptune Trojans dzhou et all d2009h ) it was shown that especially for 
larger inclinations Trojans may be in stable orbits. We have therefore undertaken an exten- 
sive numerical study of the - hypothetical - Uranus Trojan cloud for the whole range of 
inclinations 0° < i < 60°. 
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2 The method of investigation 

2. 1 The model and the grid of initial conditions 

As dynamical model for the determination of the stable regions around Uranus' L4 and 
Lg points the Outer Solar System was considered, consisting of the Sun and the four gas 
giants. The fictitious asteroids were assumed to be massless. No test computations had to 
be undertaken, because forme r studies for the Jup i ter and the Neptu n e Trojans provide d 
realistic results for that model fovorak et alJ d2008n : lzhou et al.l j2009h : lErdi et al.l d2009h l. 
Thus the equations of motion (only in the Newtonian framework) were integrated using 
the Lie-series method, a code which uses an adaptive step size and which has already bee n 
extensively tested and used in man y numerical studies (e.g. lHanslmeier and Dvorak! dl984h : 
iLichteneggerl Jl984h : lDelv«3 Jl984h l. 

Different time scales were tested: First estimations have been undertaken for only 10 6 
years of integration time to separate the regions in phase space which are more or less 
immediately unstable. Then, for the remaining stable orbits, the integrations were extended 
up to 10 7 and 10 s years, and finally for three selected orbits the time scale was set to the age 
of the Solar system (5 x 10 9 years). 

The initial conditions in the vicinity of the triangular equilibrium points were chosen 
such that the orbital elements mean anomaly, eccentricity and longitude of the node for 
the Trojans were set to the corresponding values of Uranus and only the argument of the 
perihelion was set to CO = COjj ram , s ± 60° . 

For one run different values of the semimajor axes close to the one of Uranus were taken: 

&Tto jan ~~ ®U Tonus 

± m X 0.007 AU, for 1 < m < 50. The inclination was set to be equal to 
the one of Uranus. For a second run the semimajor axis was set to the one of Uranus, but 
now the inclinations of the fictitious Trojans were changed in the range ivranus < hmjan < 
60° for steps of Ai= 1°. Both studies were undertaken separately for the leading and the 
trailing Lagrange-points (see chapter 4). For the whole grid of initial conditions around L4 
we integrated some thousand orbits to get a first picture of the stable region (see chapter 3, 
Fig.0. 

2.2 The methods of analysis 

Different tools of analysis were used to determine on one hand the stable regions in phase 
space and on the other hand to find out the reason for the instability of orbits: 

- the maximum eccentricity within the integration time, 

- the libration amplitude around the Lagrange point, 

- the escape time for the unstable orbits, 

- the involved frequencies. 

For the first three tools there is not much to say about, because they are self explaining, 
except that we define the libration width (or libration amplitude) as the difference of the 
maximum and minimum value of the mean longitude-difference of the Trojan with respect 
to Uranus. The program package SigSpec used for the frequency analysis was provided by 
Reege^ d2007l) . This is a code which is very efficient when used for data with equidistant 
time-domain sampling interval^ but also for non-equally spaced data. The computer pro- 



This does NOT mean that we used a constant step size for the numerical integration of the equations of 
motion. 
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gram is a novel method in time series analysis It incorporates the frequency and phase of 
the signal appropriately and takes also into account the prop erties of the time -domain sam- 
pling. Thus it provides a very accurate peak detection (Kallin ger et al. The benefit 
of frequency- and phase-resolved statistics and the capability of full automation are impor- 
tant also in this application with equidistant data. 

3 Stability of L\ depending on the inclination: a first overview 

To get an idea of the stability we used the results of numerically integrated orbits of several 
thousand fictitious Trojans for an integration time of 10 6 years. For an increment of 1° 
for inclinations between 0° < irmjan < 60° we used 100 orbits with different semimajor 
axes defined before. As mentioned we used three different methods for the analysis, which 
provided all more or less the same results concerning the stability of an orbit. In Fig. Q] 
we mark the amplitudes A of libration around the L4 point in different greytones and with 
contour-lines. The innermost lines on the upper part of the graph (7 > 20°) enclose orbits 
in the stable region with A < 15° (black region); the fuzzy edges show the coexistence 
of stable and unstable orbits in the transition regions. Furthermore it is interesting to note 
that for small inclinations there are no orbits with amplitude s of libration smaller than A < 
30° present, which was already mentioned in the paper by Mar zari et al ] J2003h . On the 
contrary quite small libration amplitudes occur for large inclination (42° < < 52°). The 
non existence of orbits from moderate inclinations on ( 15° < i < 3 0°) was already fou nd in 
the papers bv iNesvornvi J2000h : iNesvomv and Donesl {2002) and iMarzari et alJ J2003T) : in 
our study we can confirm it, but then, for larger eccentricities, new regions of stable motion 
appea^]. 

4 The Results of the Cuts 

The integration of thousands of orbits for longer time would need a huge amount of com- 
puter time - even with a fast program and several CPUs for running the code. We decided 
to limit the integration time to 10 7 years (or 10 8 years where appropriate) in the dynamical 
model of the OSS and specified the following initial conditions for the fictitious Trojans: the 
initial orbital elements for the massless asteroids were set to the one of Uranus, but for the 
L4 Trojans we set CO = CO + 60° and for the Lg Trojans CO = CO — 60°. To be able to deter- 
mine the sizes of stable regions with respect to the semimajor axes and the inclination of the 
Trojans we chose the following two different experiments: 

- a-cut where we varied the initial semimajor axes of 100 fictitious Trojans and set the 
initial inclination - as an illustrative example - to i = 3°; the integration time was limited 
to 10 7 years. 

- i-cut where we fixed the semimajor axis to the one of Uranus and set the inclination of 
60 fictitious Trojans to 1° < i < 60°; the integration time was set to 10 8 years. 

In addition two different runs - independently for L4 and L5 - have been undertaken to 
check any possible differences. 

2 It is the first technique relying on an analytic solution for the probability distribution produced by white 
noise in an amplitude spectrum (Discrete Fourier Transform). Returning the probability that an amplitude 
level is due to a random process, it provides an objective and unbiased estimator for the significance of a 
detected signal. 

3 compare section 4.2, Figf4] 
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contour levels: 




18.90 18.95 19.00 19.05 19.10 19.15 19.20 19.25 19.30 19.35 19.40 
initial semi-major axis [AU] 



Fig. 1 The initial semi-major axes of the fictitious object (x-axis) are plotted versus the initial inclination 
of the Trojan (y-axis). The grey-tone visualizes the libration amplitude around the Lagrange point L4; the 
contour-lines show the limits of the libration amplitudes with a step of 15° (for detail see in the text) 



4.1 The a-cuts 



As mentioned before these computation are very cumbersome from the point of view of 
the CPU times. Therefore we made the following decisions for the initial conditions of 
the a-cuts: the inclinations of the Trojans were - as mentioned before - i = 3° and the 
semimajor axes were set to 100 values between 18.9 < a < 19.6 AU, such that the Lagrange 
points should be in the middle of this interval. This choice showed the apparent 'asymmetry' 
between L4 and Lg, but displayed well the actual location of the Lagrange points. This kind 
of asymmet ry also appeared in the inves ti gations of the Jupiter Trojan s already in the first 
studies (e.g. iHolman and Wisdom! 1 1993b: iNesvornv and Donesl 1 2002)) and was discussed 
in more detail in recent pap e rs jRobutel and Gabernl d2006h : iLhotka et all ( 1200 8h ). also for 
Neptune ( Zhou etaD d2009b ; iDvorak et al.U 2007. 2008)). The apparent 'asymmetry' is due 
to the differently chosen initial conditions for the hypothetical Trojans when one investigates 
the area around L4 and Lg because of the presence of Jupiter which changes the location of 
the equilibrium points. Strictly speaking they don't even exist in the n-body problem, but 
they can be defined as the minimum in the libration angle. Nevertheless the stable regions 
are symmetric and it is sufficient to explore only one (e.g. L4). In fact comparing the stability 
diagram of the fictitious Trojans with different semimajor axes for L4 and Lg first of all the 
shift in the stable region is eye-catching: the Lagrange point itself is lying at a = 19.235 AU, 
whereas the middle of the stable area for L4 is at a = 19. 18 AU, and for Lg it is at a = 19.3 
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AU (cf. Figs.[2][3]l. The respective values for position in a for the Lagrange points were then 
chosen as the initial semimajor axis for the i-cuts described in the next section. 

The time interval to show the extent of the stable region with respect to the semimajor 
axis was set to 10 7 years; we are aware that this is not enough to determine the effective 
largeness, but gives a good second estimate after the results for the 10 6 years integration 
for the Lagrange point L4 of Fig.[TJ In fact one can see that in Fig.[2]the outer parts, which 
still show stable orbits in the 10 6 year integration, are unstable, but the central area is still 
very stable with librations with amplitudes A < 80°. It has been mentioned in the former 
section that on the edge stable and unstable orbits are close to each other which explains the 
unstable gap around a = 19.07 AU especially visible in the libration plot (Fig.[2]bottom). 

This unstable gap is not visible for orbits around the Lagrange point Lg; it shows a 
smooth stable region with increasing libration amplitudes on both sides of the equilibrium 
point itself (Fig.[3]bottom). The different escape times and larger eccentricity values in the 
unstable region are not surprising, because these orbits are - after leaving the region around 
the Lagrange point - in a chaotic state suffering from multiple encounters with the planets. 



4.2 The i-cuts 

In Fig. [TJ there are two main regions with stable orbits separated by a curved strip of insta- 
bility for an inclination around i = 16°. Although we expected that the stable regions will 
shrink for an integration time 100 times longer (10 8 years), the opening of an unstable gap 
around i = 8° was a surprise. In the next section, where we carefully discuss the dynamics 
of three fictitious Trojans, a possible explanation will be given. 

According to our results we can distinguish four different connected region in Fig. [4] 
where the orbits are stable with e max < 0.2 (upper graph), no escapes occur up to 100 million 
years (middle graph) and only with small libration angles (bottom graph): 
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Uranus Trojans L 4 , a= 19.18 AU 




10 20 30 40 50 




10 20 30 40 50 60 



initial inclination [deg] 



Fig. 4 i-cut for L4 : inclination versus the maximum eccentricity (top), versus the escape time (middle) and 
versus the libration width around the Lagrange point (bottom). 
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Uranus Trojans L 5 , a = 1 9.30 AU 




10 20 30 40 50 60 




10 20 30 40 50 60 

initial inclination [deg] 



Fig. 5 i-cut for L5: inclination versus the maximum eccentricity (top), versus the escape time (middle) and 
versus the libration width around the Lagrange point (bottom). 



A for 0° < hrojan < 7° 

Bfor9°<i rro ; fl>I <13° 
Cfor3i°< i T roj<m<36° 
D for 38° < irrojan < 50° 

This description and the finding of the regions A to D has been done for the L4 environ- 
ment. The picture for L5 is different in the following features: the unstable strip at around 
i = 8° is slightly larger, also the unstable gap close to i = 37° is broader which makes region 
D significantly smaller. What we expect is that for even longer integration time the stable 
region C will disappear and we will be left with three 'islands' of stability with respect to 
the inclinations around both Lagrangian equilibrium points. We could not identify the small 
instability strip between A and B with a special secular frequency (visible also in the signal 
of the Ura3 orbit). The instability between B, respectively C and D is caused by #5 (Jupiter) 
and g-i (Uranus) accordi ng to our analysis of the main frequencies which corresponds quite 
well with the results of iMarzari et al.l |20Q3) (their Fig. 10). No stable orbits are present 
for i > 50° which is due to the presence of gs (Neptunej^- A detailed investigation of all 
frequencies involved in these regions is in preparation, where not only cuts will be used but 
also the extensions in semimajor axes and inclination. 



4.3 Escapes during 100 Million years 

Another plot was drawn concerning the number of escapers from the region around L4 of 
Uranus during 100 million yea rs. Here we show the numbe r of escapers and compare the 
results with the ones derived in lNesvornv and Donesl ll2002l) (Fig.[6]>. In our investigations, 



4 gs to gg are the fundamental frequencies of the giant planets Jupiter (5) to Neptune (8) 
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1e+05 1e+06 1e+07 1e+08 

escape time [yrs] 

Fig. 6 Escape times of Uranus Trojans close to Z.4; the full line summarizes the results. 



where we restricted the initial semimajor axes to 19.1 < a < 19.3 AU, ~ 3100 orbits were 
integrated. Here it is possible to distinguish for different inclinations; it is evident that the 
escape behaviour is quite different for = 1,20,40 and 60°. For all the inclinations together 
the plot shows that at the beginning - the first million years - almost no objects escape, but 
then a constant depletion up to 100 million years occurs. Finally after the total integration 
ti me some 20% of the original p opulation survive, which is comparable to the results derived 
bv lNesvornv and Donesl (l2002h . The different escaping rate for different initial inclinations 
reveals that there are different mechanisms taking effect on different timescales, which will 
be analyzed in detail in our future work. 



5 Long term integrations 

To assure that the determined stable regions will survive for times up to the age of the solar 
system three orbits close to the Lagrange point L4 were integrated for 5 x 10 9 years for low 
inclined orbits. The initial conditions were the ones of Uranus except for the inclinations: 

- Ural with irrojan = tUranus, 

- Ura2 with irrojan = tUranus + 2°, 

- Ura3 with irrojan = ft/row* + 4°. 

In the respective plot where we compare the eccentricities of fictitious Trojans directly 
with Uranus (Fig. [7}, it is evident that there is not a big difference in the time development 
of the eccentricities; only the most inclined Trojan Ura3 suffers from larger amplitudes 
than Ural and Ura2 in e. The comparison of the inclination plots unveils why there is an 
instability gap at about i = 6°: Whereas on the first graph (Fig. [Sj> no difference is visible 
at all between the Trojan Ural and Uranus itself, already in the graph for Ura2 (Fig.|9]l a 
significant feature can be observed. Large periodic changes occur with i max > 4°. From the 
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Fig. 7 The eccentricities of 3 fictitious Uranus Trojans after 5 billion years compared to the one of Uranus: 
Ural (top), Ura2 (middle) and Ura3 (bottom). 



Ura 1 for first 50 million years 




4.95e+06 4.9Ge+06 4.97e+GG 4.98e+0S 



Fig. 8 Inclination of Uranus and Ural for the first 50 and the last 50 million years of integration (upper, 
respectively lower graph). Note that the signal for Ural and Uranus are undistinguishable. 
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Ura 2 for first 50 million years 




Fig. 9 Inclination of Uranus and Ura2 for the first 50 and the last 50 million years of integration. 
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Ura 3 for last 50 million years 



: l_ i i i 

4.95e+06 4.96e+G6 4.97e+0G 4.93e+06 4.99e+G6 5e+ 




Fig. 10 Inclination of Uranus and Ura3 for the first 50 and the last 50 million years of integration. 
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Fig. 11 The two main periods in the element k = ecos(I2 + o) for the 5 billion years integration of Ura3 
namely the 4.4 x 10 s years period and the period which changes between ~ 10 7 and 2 x 10 7 years compared 
to the 'same' period of ~ 10 7 of Ura2 which is shown as as straight WALL in the middle of the graph: periods 
(y-axis) are plotted versus the time (x-axis) and the amplitudes (z-axis). 



graph of Ura3 (Fig.llOt one can see that the period of the variation of the inclination at the 
end of the 5 billion years is only half of the period at the beginning. 

A carefully undertaken frequency analysis showed that in contrast to Ural and Ura2 for 
Ura3 a dramatic change in the main frequencies occured during the dynamical evol ution. 
This i s a sign of ch aotic behaviour as was already pointed out in many papers (e.gjEaskar 
dl990l . ll994lll997t) ). To visualize this change we used the technique of running windows, 
each section spanning 50 million years with an overlap of 25 million years. The changes 
are well visible for the larger periods of some 12 to 18 million years which also show 
significantly larger amplitudes than the second largest period of about 4.4 x 10 years (Fig. 
CCD. 

In addition, at the lower end of the high inclination window, we made numerical inte- 
grations up to lGyr which showed that even orbits close to the unstability window have the 
typical behaviour of changing the periodicity. This can be seen already after several hundred 
million years in Fig[l2]for two orbits, one is well inside the stability region (z = 42°) and 
the other one already close to the edge (i = 40°). In fact the last one becomes unstable after 
another 100 Myrs (at 946 Myrs); this behaviour gives notice already when one compares the 
dynamics for the first 200 Myrs and the last 200 Myrs of this hypothetical Trojan (compare 
the 2 nd and 4 th panel in Fig. [Till. 



6 Conclusions 

In our study we explored the regions of possible Trojan asteroids of Uranus for a large 
range of their inclinations up to = 60°, which has never been done before. The method 
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Fig. 12 Eccentricity of two orbits close to the lower instability border for inclinations / = 42° (I s ' and 3 
panel) and for i = 40° (2 nd and 4"' panel) for the first 200 Myrs (upper panels) and from 700 Myrs to 900 
Myrs (lower panels); time is measured in Myrs. 



of determining these regions was the analysis of the results of straightforward numerical 
integrations of the full equations of motion in the model of the OSS. Besides the escape from 
the region we also checked the amplitudes of libration of these fictitious asteroids around 
the Lagrangian equilibrium points. To assure the survival of the 'stable' Trojans we did 
also integrations of three specially chosen Trojans on low inclined orbits for the age of the 
Solar system of 4.5 x 10 9 years, where the accurate analysis of the frequencies involved was 
undertaken. Our results showed that close to the first instability window 6° < irrojan < 8° 
the strong variation of the frequency of a special secular resonance indicates chaotic motion. 
Two other unstable windows with respect to the inclination are 14° < irrojan < 31° and 
around irrojan = 37°. From irrojan > 50° (for L4) and irrojan > 45° (L5) no stable Trojan 
orbits exist for Uranus. These stable islands in inclination are surrounded by more or less 
small stable regions in semimajor axis. 

According to the results of iMorbidelli et afl d2005h Jupiter Trojans could be trapped 
during the phase of early migration of the planets in our Solar system. If that was also the 
case for Neptun^l and Uranus Trojans (and why not) we should be able to detect them. So 
the question is: why don't we observe them? 
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